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AT-Jettiness: An Inclusive Event Shape to Veto Jets 
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Jet vetoes are essential in many Higgs and new-physics analyses at the LHC and Tevatron. The 
signals are typically characterized by a specific number of hard jets, leptons, or photons, while the 
backgrounds often have additional jets. In such cases vetoing undesired additional jets is an effective 
way to discriminate signals and background. Given an inclusive event sample with A*' or more jets, 
the veto to have only A*' energetic jets defines an "exclusive" A^'-jet cross section. This strongly 
restricts the phase space of the underlying inclusive A^'-jet cross section and causes large double 
logarithms in perturbation theory that must be summed to obtain theory predictions. Jet vetoes 
are typically implemented using jet algorithms. This yields complicated phase-space restrictions and 
one often relies on parton-shower Monte Carlos, which are limited to leading-logarithmic accuracy. 
We introduce a global event shape 'W-jettiness" , rjv, which is defined for events with N signal jets 
and vanishes in the limit of exactly A' infinitely narrow jets. Requiring rjv ^ 1 constrains radiation 
between the A'' signal jets and vetoes additional undesired jets. This provides an inclusive method 
to veto jets and to define an exclusive A'^-jet cross section that can be well-controlled theoretically. 
A'^-jettiness yields a factorization formula with inclusive jet and beam functions. 



Introduction. At the LHC or Tevatron, hard interac- 
tions involving Higgs or new-physics particles are identi- 
fied by looking for signals with a characteristic number of 
energetic jets, leptons, or photons The backgrounds 
come from Standard Model processes producing the same 
signature of hard objects possibly with additional jets. 
An example are top quarks decaying into W plus 6-jet, 
which is a major background for H — )■ WW [2[. When 
reconstructing masses and decay chains of new-physics 
particles additional jets can cause large combinatorial 
backgrounds. Standard Model processes can also fake 
a signal when a jet is misidentified as lepton or photon, 
a typical example being H — > 77. 

Thus, a veto on additional undesired jets is an effective 
and sometimes necessary method to clean up the events 
and discriminate signal and the various backgrounds. 
More generally, one would like to measure an "exclusive" 
A^-jet cross section, pp — >■ XL{Nj), to produce N signal 
jets j where the remaining X contains no hard (central) 
jets. Here, > and L denotes the hard leptons or 
photons required as part of the signal. 

We introduce an inclusive event shape "A^-jettiness" , 
denoted r^r and defined below in Eq. ([T]). For an event 
with at least N energetic jets, tn provides an inclusive 
measure of how A^-jet-like the event looks. In the limit 
Tat — > the event contains exactly N infinitely narrow 
jets. For tn ^ 1 the event has hard radiation between 
the N signal jets. Requiring tat <C 1 constrains the ra- 
diation outside the signal and beam jets, providing an 
inclusive way to veto additional central jets. It yields 
an inclusive definition of an exclusive A^-jet cross section 
with a smooth transition between the case of no jet veto, 
tn ^ 1, and the extremely exclusive case tat — > 0. 

Vetoing additional jets imposes a phase-space restric- 
tion on the underlying inclusive A^-jet cross section to 
produce A^ or more jets with the same L. Irrespective of 
its precise definition, the jet veto introduces a jet resolu- 



tion scale fij that characterizes this restriction, i.e. the 
distinction between A^ and A^-l-l jets. Hence, the exclu- 
sive A^-jet cross section contains phase-space logarithms 
a" In™' {iij / fijj) , where m < 2n and /i// is the scale of the 
hard interaction. For r^r, ^'j/^'j^ ~ tjv ^ 1. Generically 
there is always a hierarchy fij fj,H, which becomes 
larger the stronger the restrictions are. These large log- 
arithms must be summed to obtain reliable predictions. 

Jet vetoes are typically implemented by using a jet al- 
gorithm to find all jets in the event and veto events with 
too many energetic jets. Jet algorithms are good tools 
to identify the signal jets. However, they are not nec- 
essarily well-suited to veto unwanted jets, because the 
corresponding phase-space restrictions are complicated 
and depend in detail on the algorithm. This makes it 
difficult to incorporate the jet veto into explicit theoret- 
ical calculations and inhibits a systematic summation of 
the resulting large logarithms. In this case, usually the 
only way to predict the corresponding exclusive A^-jet 
cross section is to rely on parton shower Monte Carlos 
to sum the leading logarithms (LL). For particular jet 
algorithms, the resolution 1/23 defines the transition from 
2 to 3 jets. Next-to-leading logarithms for this and other 
hadron-collider event shapes were summed in Ref. (3|- 

Vetoing jets by cutting on an inclusive variable like tat 
has several advantages. First, we can go beyond LL or- 
der, because the logarithms from the phase-space restric- 
tion, a" In^TAT, are simple enough to allow their system- 
atic summation to higher orders. Moreover, the theory 
predictions with factorization can be directly compared 
with experiment without having to utilize Monte Carlos 
for parton showering or hadronization. Experimentally, 
Tat reduces the dependence on jet algorithms and might 
help improve the background rejection. 

Definition. A^-jettiness is defined as 
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FIG. 1: Dilferent situations for the application of N-jettiness 



Jet 2 2 / Jets 3 a 
(c) pp — > leptons plus jets. 



As we discuss below, this definition of tn yields a fac- 
torization formula with inclusive jet and beam functions 
and allows the summation of logarithms to next-to-next- 
to-leading logarithmic (NNLL) order. The sum over k 
in Eq. ([T]) runs over the momenta pk of all measured 
(pseudo-)particles in the final state excluding the signal 
leptons or photons in L. (Any other leptons or photons, 
e.g. from hadronic decays, are included in the sum.) For 
simphcity we take all pk to be massless. The qa, qt, and 
qi, qN are a fixed set of massless reference momenta 
for the two beams and the N signal jets, 
1 



(1, 



q^j^Ej{l,nj), J = {l,...,N}. (2) 

The Ej and hj correspond to the energies and directions 
of the N signal jets (for both massive and massless jets). 
Their choice is discussed below. The beam reference mo- 
menta qa and qb are the large momentum components of 
the colliding partons along the beam axis (taken to be 
the z axis). They are defined by 
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• (gi H h qAT + q) , 



(3) 



and analogously for Xh with a ^ b. Here, q is the to- 
tal momentum of the non- hadronic signal L. In Eq. ([T]), 
= XaXbE^j^ is the hard interaction scale, and the dis- 
tance of a particle with momentum pk from the jets or 
beams is measured by q„i ■ p^. If L contains missing en- 
ergy, so q and Xa^b are not known, one can use a modified 
distance measure as we discuss below Eq. pD) . 

The minimum for each k in Eq. ([T]) associates the par- 
ticle with the closest beam or jet, appropriately dividing 
the hadronic initial-state radiation (ISR) and final-state 
radiation (FSR). Soft particles and energetic particles 
near any jet or beam only give small contributions to the 
sum. For 2 — >• TV scattering of massless partons, tn — 0. 
Energetic particles far away from all jets and beams give 
large contributions. Hence, for rjv ^ 1 the final state has 
N jets, two forward beam jets, and only soft radiation 
between them. In this limit Xa,b are the momentum frac- 
tions of the annihilated partons, and Y = \n{xa/ Xb) /2 is 
the boost of the partonic center-of-mass frame. 



N — 2 for e^e^ — >■ jets. In e~^e~ collisions there is no 
hadronic ISR, so we drop the qa,b ■ Pk entries in Eq. ([1]). 
Now is the total invariant mass of the leptons and Y = 
0. In the two-jet limit, the jet directions are close to the 
thrust axis i, defined by the thrust T ~ max^ \i-pi\/Q. 
Hence we can choose 



q^^lQ{l,i). 



q^^lQ{l,-i) 



(4) 



as reference momenta, and Eq. ([T]) becomes 
1 
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Eknim{l — cosOk, 1 + cosOk} , (5) 



where 9k is the angle between pk and t. The minimum 
divides all particles into the two hemispheres perpendic- 



ular to i as shown in Fig. 1(a) For r|*^ <^ 1, the total 



invariant mass in each hemisphere is much smaller than 
Q, so the final state contains two narrow jets. In this 
limit, tI'^ — 1 — T, and a factorization theorem exists for 
dcr / drf"^ , which can be used to sum logarithms of t|'^ 0] . 
For a given jet algorithm with resolution parameter y, 
the value 1/23 marks the transition between 2 and 3 jets. 
Thus requiring 1/23 <C 1 also vetoes events with > 2 jets. 

iV = for Drell-Yan. Next, consider the isolated 
Drell-Yan process, pp — X£'^i~ with no hard central 



jets, shown in Fig. 1(b) We now have ISR from the in- 



coming partons, but no FSR from jets. From Eq. Q we 
have 
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where q^ and qT are the dilepton invariant mass and 
transverse momentum, and Y equals the dilepton rapid- 
ity. Now, Q"^ — q^ -\- q^ and Eq. ([T]) becomes 



feT|min{e^~''^ 6-'^+'"=}. 



(7) 



where \pkT\ and r]k are the transverse momentum and 
rapidity oi pk- The qa and qb dependence in Eq. ([T]) ex- 
plicitly accounts for the boost of the partonic center-of- 
mass frame. For y = 0, the minimum in Eq. ([7]) divides 
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all particles into two hemispheres perpendicular to the 
beam axis (analogous to t above). For Y 0, the hemi- 
spheres are boosted with their boundary now at Y, and 
the beam jet in the direction of the boost is narrower 
than the other, as depicted in Fig. 1(b) Contributions 



to To from large rapidities are exponentially suppressed 
by |p/cT|e~''"'' « 2iJfce~^l'"=l, so particles beyond the de- 
tector's rapidity reach give negligible contributions. 

Beam thrust [1] is given by tb = \/T 
obtained by choosing Xa,bEcni - 



q^/q^TQ. It is 



q2g±Y case q^ and 
Y are measured rather than the longitudinal components 
Ua^b ■ q in Eq. For tq <C 1 the hadronic final state 
can only contain soft radiation plus energetic radiation 
in the forward directions, so {qxl ^ Q and tb = tq. A 
factorization theorem for dcr/drs at small tb was derived 
and used to sum logarithms of tb to NNLL in Refs. 

General case. For pp — ^ XL{Nj) we have both ISR 
and FSR. We select candidate signal events by measuring 
L and running a jet algorithm to find the N signal jets 
and their momenta pj. The conditions on the jets and L 
that define the signal are encoded in the cross section by 
a measurement function Fn{{pj}, L). Generically, Fjv 
will enforce that there are at least N energetic jets that 
are sufficiently separated from each other and the beams. 
We now use the measured jet energies and directions to 
define the massless reference momenta qj in Eq. 
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nj =pj/\pj\ 



(8) 



while qa and qt are given by Eqs. ^ and ([3]). 

Taking the minimum in Eq. ([IJ combines the previous 
cases in Eqs. ([5]) and ([7]). It divides all particles into 
jet and beam regions that are unique for a given set of 
reference momenta and whose union covers all of phase 



space, as illustrated in Fig. 1(c) The boundary between 



any two neighboring regions is part of a cone and is such 
that the sum of the total invariant masses in the regions 
is minimized (or in case of a beam region the virtuality 
of the incoming colliding parton). 

For events with small r^r all jet algorithms should agree 
how energetic radiation is split up between the jets and 
beams, and only differ in their treatment of softer par- 
ticles. Thus, they all give the same rij and Ej up to 
power corrections, while the split up of the soft radia- 
tion is determined by tn itself. Hence, the dependence 
of tn on the jet algorithm is formally power suppressed. 



alg.l 



alg.2 



T^'° = Tjy"' + 0{t^), as seen in Eq. (jH)) below. 



To measure r^v, we still rely on having a suitable jet 
algorithm to find the N signal jets but not more so than 
if we were not measuring Tpf. Imagine the jet size in the 
algorithm is chosen too small such that the algorithm 
divides what should be a single signal jet into several 
narrow jets ll[. In this case, the jet algorithm yields a 



poorly reconstructed signal irrespective of measuring tn- 
Since the jet veto is now provided by tat, this situation 
can be avoided because we do not have to rely on the jet 



algorithm to identify additional jets and so can use an 
algorithm that can be forced to always yield at most N 
jets. This is in fact the most natural thing to do when 
one is looking for N jets. Therefore, using tjv as jet veto 
could also help improve the signal reconstruction. 
Generalizations. We can generalize tn to 

= '^mm{da{pk),db{pk),di{pk), . . . ,dN{pk)} , (9) 



where dm{pk) can be any infrared-safe distance measure. 
In Eq. (H]), drn{pk) = Sq^ ■Pk/Q'^ with 

2qj -pk = \pkT\ \qjT\ (2coshA77,/fc - 2cosA(/),/fc) . (10) 

Here, Arjjk and A(f)jk are the rapidity and azimuthal 
distances between qj and pk. If these are small, the factor 
in brackets reduces to the familiar = (A?])^ + (A(/))^. 

Different measures that are boost-invariant along the 
beam axis can be obtained by modifying the dependence 
on rapidity, \qjT\, and Q in Eq. (fTO|) . A geometric mea- 
sure, which is independent of |(j!/t|, is 



db{pk) 



_ \PkT[ ^^Y+v, 



dj{pk 



\PkT\ 



(2 cosh Aryjfe - 2 cos A^jfc) . 



(11) 



It evenly divides the area rather than invariant mass be- 
tween neighboring regions, such that more energetic jets 
also get more invariant mass. 

If L contains missing energy, then Xa,b in Eq. ([3]) and 
thus Q and Y are not known. For Q, one can use any 
hard scale, like the \qjT\ of the hardest jet or leave it out, 
since it only serves as an overall normalization. In the 
beam measures da^b{Pk) we can simply set y = 0, which 
defines them in the hadronic center-of-mass frame. 

A''-jettiness does not split events into N , + 1, A^ + 2, 
etc. jets like a traditional jet algorithm. But we can con- 
sider using Tjv to define an "exclusive A^-jet algorithm" as 
follows: First, we use a geometric measure and find the 
directions hj and boost Y that minimize r^r, analogous 
to finding t for e+e" — >■ jets. This might actually allow 
one to get an estimate of Y even in the case of missing 
energy by exploiting the asymmetry in the beam jets. 
Second, we determine the jet energies by summing over 
the particles in each jet region. (To reduce the sensitivity 
to the underlying event and pile-up, one can weigh the 
sum over energies by the distance from fij.) 

Factorization formula. We now use t^t again as de- 
fined in Eq. ([J). For tat < 1, QCD ISR and FSR can be 
described in soft-collinear effective theory [6] at leading 
power by A^ + 2 independent sectors for coUinear parti- 
cles close to each q^ with m = {a, 6, J} and a separate 
sector for soft particles. By power counting, J-coUinear 
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particles are closest to qj, so for the J-collinear sector 
m\n{2qm ■ Pk) = ^ 2qj ■ pk = sj , (12) 

where (up to power corrections) sj is the total invariant 
mass in the J-collinear sector. Similarly, the sum over 
the beam collinear sectors yields the total (transverse) 
virtuality of the colliding partons, ta and tb. Therefore, 



ta + h + '^sj + ^ mm{2q,n ■ Pk} ■ (13) 



A: soft 



The sum in the last term is now restricted to the soft 
sector. Combining Eq. (|T3l) with the analyses in Refs. 0, 
0| yields the factorization formula for A^-jettiness [l^l 

^= JdXadXb jdSd^Liq) Jd^N{{qj}) 

X FNi{qm}, L) {2^f5''[qa + qb-Y.qj- q) 

J 

X Hij^f,{{qm},L,^i)Y\_ J dsj Jf,j{sj,ii) 

X dtaBi{ta,Xa,fJ.) dtb Bj{tb, Xt, fJ.) 



xSr-{rN-'-^^^^'^fi^,Um},t^). (14) 



iij^K({q,„}, L) contains the underlying hard in- 



Here, Hi. 

teraction i{qa)j{qb) L{q)Ki{qi) ■ ■ ■ KpfiqN), where i, j, 
and Kj denote parton types, and the sum over ij, k is 
over all relevant partonic channels. It is a matrix in color 
space given by the IR- finite parts (in pure dim. reg.) of 
the squared partonic matrix elements in each channel. 
The iV-body phase space for the massless momenta qj is 
denoted d$„({gj}), and that for L by d^L{q)- 

The inclusive jet and beam functions, Jkj{sj) and 
Bij{ta,b,Xa,b), describe the final and initial state radia- 
tion emitted by the outgoing and incoming partons from 
the hard interaction. The latter also determine the mo- 
mentum fractions Xa,b of the colliding partons and are 
given by parton distribution functions ff (^, /i) as [^, 0] 



B,{t,x,fi) = ^ 



■{t,^,f?jU{C,fi). (15) 



The lii' are perturbative coefficients describing collinear 
ISR, and at tree level Bi{t, x, ^) = S{t)fi{x, fi). The last 
term in Eq. (jl3p is the contribution to tjv from soft par- 
ticles in the underlying event. It is described by the soft 
*'*(t^^*, {<Zm}), which depends on the jet's 



function 



angles hi ■ hm and energy fractions Ei/ Em- Like H, it is 
a color matrix, and the trace in Eq. (fT4|) is over tr(iJS'). 

In Eq. (fT4|). all functions are evaluated at the same 
renormalization scale /i. Large logarithms of r^r in 
da/dTM are summed by first computing H{ph), Jip-j)! 



B{fj.B), S{fis) at the scales fin — Q, IJ-j — fJ-B — y/mQ, 
fis — tnQ, where the functions contain no large loga- 
rithms, and then evolving them to the scale /i. This evo- 
lution is known analytically Q and the required anoma- 
lous dimensions are aheady known to NNLL be- 
cause we have inclusive jet and beam functions. NNLL 
also requires the 0{as) corrections for each function, 
which are known for J and B. The 0{as) hard function 
is determined by the one-loop QCD matrix elements. For 
tn ^ Aqcd/Q, S^is) can be computed perturbatively 
and will be given in a future publication. 
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This can be tested by comparing the total energy in each 
region defined by rjv with the energy from the jet algo- 
rithm. If these are very different, but at the same time rjv 
is small, then there are additional energetic particles near 
the signal jets that the algorithm should have included. 
Here, _Fjv enforces distinct collinear sectors with 1 — n; ■ 
fim. ^ Tjv and Em/Q ^ TN. We assume _Fjv only de- 
pends on the large components q,j of the jet momenta, 
Pj = qj[l -I- O(rjv)], and that L only couples to the QCD 
subprocess via a hard interaction. We also assume that 
Glauber gluons do not spoil this factorization. 



